Towards a New Proof of Anderson Localization 
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The wave function of a non-relativistic particle in a periodic potential admits oscillatory solutions, 
the Bloch waves. In the presence of a random noise contribution to the potential the wave function is 
localized. We outline a new proof of this Anderson localization phenomenon in one spatial dimension, 
extending the classical result to the case of a periodic background potential. The proof makes use of 
techniques previously developed to study the effects of noise on reheating in inflationary cosmology, 
employing methods of random matrix theory. 



I. INTRODUCTION 

A classic problem in non-relativistic quantum mechan- 
ics is the propagation of a particle, e.g. an electron, in a 
periodic potential, set up e.g. by a lattice of ions. It is 
known that the time-independent Schrodinger equation 
for this problem admits modulated periodic solutions, the 
so-called Bloch waves In the real world, however, 
we expect that the potential will not be perfectly peri- 
odic due to the presence of various types of disorder, e.g. 
thermal noise. As was first shown by Anderson 0, if the 
noise is modeled as a random contribution to the poten- 
tial, then for sufficiently large disorder the wave functions 
become localized. In one spatial dimension, a stronger re- 
sult holds: it can be shown J3] that for any amount of 
disorder, the wave functions become exponentially local- 
ized, i.e. instead of being oscillatory, the solutions now 
are exponentially decaying in space about the location of 
the center of the wave packet. This result was extended 
to two spatial dimensions in . Anderson localization is 
reviewed in 0, 0, 0] ■ A mathematically rigorous proof of 
one dimensional Anderson localization for constant back- 
ground potential is given in [1| , and there is an extensive 
mathematical literature on the subject. 

As is well known, in one spatial dimension the time- 
independent Schrodinger equation can be mapped into 
a classical time-dependent wave equation by a simple 
variable transformation in which the wave function ^ 
becomes a classical field x ^nd the spatial variable x 
is transformed to time t. After this mapping, the 
Schrodinger equation for the wave function ^' in a po- 
tential which is the sum of a periodic term Vp{x) and a 
random noise term Vr{x) becomes the relativistic wave 
equation (in momentum space) for a scalar field with a 
mass which contains one contribution to the mass which 
is periodically varying in time, and a second contribution 
which is random in time. 

The inflationary scenario [§] (see also for a recent 
review emphasizing both successes and problems of infla- 
tion), the current paradigm of early universe cosmology. 



is based on the dynamics of classical scalar fields coupled 
as matter source to Einstein's theory of General Relativ- 
ity. According to the inflationary universe scenario, there 
is a period in the very early universe in which space ex- 
pands almost exponentially. This is obtained by making 
use of a scalar field Lp whose energy functional is domi- 
nated by a potential energy density contribution which is 
almost constant in time. Regular matter can be modeled 
as a second scalar field x- 

During the period of inflation, an exponentially in- 
creasing fraction of the energy density of the universe is 
stored in the held Lp. Hence, to make contact with what is 
observed today, a phase at the end of the period of infla- 
tion when the energy transfers from ip to regular matter 
is crucial. This phase is called the "reheating stage" , and 
works in the following way. The period of inflation ter- 
minates once Lp begins to oscillate about the minimum of 
its potential. Due to a resonant coupling between ip and 
the matter fleld x, energy can be transferred from ip to 
X- The theory of reheating is obtained by studying the 
evolution of x in the presence of an oscillating (p field. 

The equation of motion for x is the Klein- Gordon equa- 
tion with a correction term due to the coupling with (p. 
Treating x at the linearized level, each Fourier mode of x 
evolved independently. Neglecting the effects of the ex- 
pansion of the universe \Tf\ , and for the specific coupling 
between the two fields which we specify below, the equa- 
tion of motion for such a Fourier mode of x is that of a 
harmonic oscillator whose mass has two time-dependent 
contributions, one periodic in time (due to the oscillatory 
dynamics of Lp during the reheating period, the other be- 
ing an aperiodic random noise term (due e.g. to quantum 
fluctuations in Lp). 

In the absence of noise, the equation of motion is the 
Mathieu equation and falls into the class o f eq uations 
first studied by Hill, Floquet and Legendre [ll|. As a 
function of the value of the bare mass (the term in the 
mass independent of time), there are stability bands (for 
which the solutions are oscillatory) and instability bands 
(for which the solutions are characterized by overall ex- 
ponential growth or decay. The coefficient describing the 
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exponential growth is called the "Lyapunov exponent" . 
Applied to the case of cosmology, then in the absence 
of noise there are bands of Fourier modes for which the 
mass lies in the instability band and for which there is 
resonant increase in the amplitude of x, which in terms of 
phy sics corresponds to resonant production of particles 

Not too long ago pj|, the effects of a particular type 
of random noise on the resonant production of particles 
during inflationary reheating was studied. It was shown 
that the Lyapunov exponent /ife(g) which describes the 
exponential growth of each Fourier mode k of the scalar 
field in the presence of the noise q is strictly larger than 
the corresponding Lyapunov exponent ^fe(O) in the ab- 
sence of noise |28j . In particular, this implies that in 
the presence of noise, every Fourier mode grows expo- 
nentially. The stability bands which are present in the 
absence of noise disappear. 

In this Letter, we review the above analysis, trans- 
late the results into the language of the time- independent 
Schrodinger equation, and in this way immediately ob- 
tain a new proof of Anderson localization which extends 
to the case of a periodic background potential. 

II. REVIEW OF RESULTS ON THE EFFECTS 
OF NOISE ON REHEATING IN INFLATIONARY 
COSMOLOGY 

We begin by reviewing the results of [3] ■ We consider 
the simplest model which describes both the inflation- 
ary phase of the very early Universe and the period of 
"reheating" which terminates the phase of inflation and 
during which the energy of the Universe is transferred to 
regular matter. This model contains two scalar matter 
fields, if and x- The first, if, has a large potential energy 
which is changing only very slowly in time, dominates 
the total energy density of the Universe and thus gives 
rise to inflation 9]. The second scalar field, x, represents 
regular matter. It is in its vacuum state initially and 
gets excited by a coupling with ip during the final stages 
of inflation when tp(t) oscillates about the minimum of 
its potential [29]. It is of great interest in cosmology to 
study the dynamics of matter production at the end of 
inflation. 

We will assume that the interaction Lagrangian density 
takes the form [sO] 

Ant = ^gvx'^ , (1) 

and that the Lagrangian for x alone is that of a free 
scalar field with mass m^. We will also consider (p to 
be spatially homogeneous. In the absence of quantum 
fluctuations, this is a reasonable assumption since during 
the period of inflation, spatial fluctuations are red-shifted 
exponentially. In this case, the equation of motion for x 
can be solved independently for each Fourier mode. In 
the absence of expansion of the Universe, the resulting 



equation is [3l[ 

Xk + H + 9V=it)] Xk = 0, (2) 

where uj'^ = k'^ + m^. 

At the end of the period of inflation, the scalar field Lp is 
oscillating about its ground state value, which we without 
loss of generality can take to be (/9 = 0. The frequency 
u) of the oscillation is set by the mass of ip. In this case, 
the equation ^ has the form of a Mathieu equation. As 
is well known (see e.g. [Tl|, [13, [2l|), there are bands 
of values of Wk for which Xk grows exponentially. The 
exponential growth is governed by a Lyapunov exponent 
which can be extracted from the time evolution of 
Xk{i) as follows: 

/Xfc = \imt^oo^log\xk{t)\ ■ (3) 

The time scale of exponential growth is in many con- 
crete inflationary models much shorter than the typical 
expansion time of the Universe, thus justifying the ne- 
glect of the expansion. Thus H^, the parametric insta- 
bility leads to a very cfhcient energy transfer [32.] 

However, in cosmology we expect fluctuations of ther- 
mal or quantum nature to be super-imposed on the ho- 
mogeneous oscillation of (p. In fact, we believe that quan- 
tum vacuum fluctuations of ^p in inflationary cosmology 
are the seeds for the inhomogeneities in the galaxy dis- 
tribution and the anisotropies in the temperature of the 
cosmic microwave background observed today. Thus, it 
becomes important to study the sensitivity of the para- 
metric resonance instability of ([2]) in the presence of an 
oscillating Lp{t) to the presence of random noise in 'p{t). 
In [13] , this problem was studied for homogeneous, ape- 
riodic noise, modeled as the addition of a stochastic con- 
tribution to the time dependent mass. More precisely, 
we studied the equation 

Xk + K + + Xfe = , (4) 

where p is a periodic function with period 27r and q{t) is 
an aperiodic, random noise contribution [33| 

To derive rigorous inequalities on the magnitudes of 
the Lyapunov exponent with and without noise, it is nec- 
essary to make certain assumptions on the noise q{t). 
We consider noise drawn from some sample space $7 = 
C{TV) with a translationally invariant probability mea- 
sure dP{K), where k labels an element in Q, and assume: 

• The noise is ergodic, i.e. the ensemble average of 
a function of the noise equals the time average for 
almost all realizations of the noise. 

• The noise is uncorrelated in time on scales larger 
than r, the period of the oscillatory contribution 
to the mass in 

• Restricting the noise to the time interval < t < T, 
the samples q{t; k) : < t < T within the support 
of the probability measure fill a neighborhood of 
the origin in Q. 
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Given these assumptions, we were able to use a theorem 
of random matrix theory (apphed to the transfer matrix 
corresponding to the differential equation ([U) to prove 
that for almost all realizations of the noise 



^l{q) > m(0) 



(5) 



for all values of the momentum k. In particular, this im- 
plies that the stability bands of the Mathieu equation dis- 
appear once noise of the type considered here is added to 
the system. The noise in fact strengthens the instability. 
Note that the result ([5]) was proven non-perturbatively 
in [l^ for values of k in the stability bands, the case 
of interest to us here. We conjecture that the result is 
true in general, and have been able to show this at least 
perturbatively p3 |. 



III. APPLICATION TO ANDERSON 
LOCALIZATION 

Let us start from Eq. (|4]) and transform variables. 
We replace the time coordinate i by a spatial coordi- 
nate X (space being infinite, i.e. TV), and substitute the 
field variable x by the variable ip representing a wave 
function. With these substitutions, Eq. (|4]) becomes 
the time-independent Schrodinger equation for the wave 
function ip in one spatial dimension for a system with a 
potential which is the superposition of a periodic piece 
Vp{x) (coming from p{ujt) in ([4]) and a random noise piece 
Vii{x) coming from q{t) in 



Hijj = E%l] 



with 



H = - — — + V,{x) + Vn{x) 



(6) 



(7) 



The energy eigenvalue E is given by the momentum k of 
(HI). We use units in which h — 1. 

In the absence of noise, the spectrum of the 
Schrodinger equation ^ consists of bands of continuous 
spectrum, for which bounded quasi-periodic (modulated 
periodic) solutions exist, alternating with instability gaps 
(resonance intervals) in which the solution behavior is 
exponential. In both regions the solutions, the so-called 
Bloch waves [l|, exhibit Floquet behavior, namely they 
are of the form ip{x) = exp(±(/j,B -|- iaE)x)p{x), where 
p{x) is periodic of the same period as the potential Vp. 
The Floquet exponent itie = i^J■E + icts) has real part 
the Lyapunov exponent and imaginary part the rotation 
number of the solution ip{x). The Lyapunov exponent 
vanishes in the stable bands, while the rotation number 
takes (a fixed multiple of) integer values in the gaps. In 
the classical field theory problem these bands of spec- 
trum correspond precisely to the stability bands of the 
Mathieu equation (g]) (for q{t) = 0). 

Let us now turn on noise satisfying the conditions listed 
in the previous section and look for solutions of (O. Due 



to the translation invariance of the problem, we can take 
X = to be the point at which -0 takes on its maximum. 
In the language of the quantum mechanics problem, the 
Lyapunov exponent can be extracted as follows: 



lim 



1 



-log\i}}{x)\ 

X 



(8) 



Since the wave function in quantum mechanics must be 
normalized, no exponentially growing solutions are al- 
lowed. An eigenstate corresponds to a wave function 
which does not grow either for x — > oo or a; — )■ — oo. 
Thus, a positive real Lyapunov exponent corresponds 
to an exponentially decaying wave function, whereas an 
imaginary Floquet exponent with absolute value 1 corre- 
sponds to Bloch waves. A standard mathematical argu- 
ment |§, ll^j [Hi HEI shows that the eigenstates are dense 
and the spectrum is pure point. 

If we consider an energy eigenvalue E which lies in a 
conduction band (i.e. an energy band for which Bloch 
wave solutions exist) in the absence of noise, then our 
result ([5]) implies that as soon as a random noise term 
is added to the potential in ([7]), the Lyapunov expo- 
nent becomes positive, and that thus the corresponding 
wave function are exponentially decaying. This means 
that the addition of any noise satisfying the conditions 
listed in the previous section exponentially localizes the 
wave function. The localization length is inversely pro- 
portional to the Lyapunov exponent. 



IV. DISCUSSION 

We have shown that the time-independent Schrodinger 
equation for the non-relativistic Hamiltonian ([7]) does 
not admit Bloch wave solutions for any amplitude of the 
noise. Given a wave function centered at some point 
X (which we can without loss of generality take to be 
X = 0), we have shown that the wave function decays ex- 
ponentially away from x = Q. Thus, we have given a new 
proof of the classic result [s*] of exponential localization of 
states in one spatial dimension which extends the known 
result to the case of a periodic background potential. 

Our methods are based on translating the time- 
independent quantum mechanical problem to a classical 
dynamical systems problem, and thus are only applica- 
ble to the case of one spatial dimension and cannot be 
used to study Anderson localization in a higher number 
of spatial dimensions. 

The dynamical systems problem turns out to be a 
problem concerning the effects of noise on the resonant 
production of particles during the reheating phase of in- 
flationary cosmology. In mathematical language, our 
problem concerns adding aperiodic noise to the mass 
term in the Mathieu equation. We have shown that the 
Lyapunov exponent increases for any amplitude of the 
noise. Thus, the solutions exhibit exponential behavior 
both forwards and backwards in time (starting from the 
initial time t = corresponding to x = for the quantum 
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problem). Normalizability of the wave function in the 
quantum problem implies that the eigenstates are pre- 
cisely those solutions which decay away from the source 
point. 
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This is self-consistent if the period of energy transfer 
turns out to be short compared to the time scale of the 
expansion of space. 

Non-perturbatively, this was only shown for values of k 
in the stability bands. Perturbatively, the statement also 
holds for values of k in the instability bands What 
is relevant for this paper is the result for values of k in 
the stability bands. 

This toy model was already used in the flrst studies of 
reheating in inflationary cosmology [IE, EBj- 
The quantity (? is a constant which has dimensions of 
mass. It is the coupling constant. 

The exp ansion of the Universe can be taken into account 
[itI . llSl . [l9l | without affecting the result that there are 
exponential instabilities. 

At some point, the back- reaction of the energy in the x 
field on the dynamics of space-time will become impor- 
tant, and this could terminate the energy transfer [22[. 
The mathematically much more complicated problem of 
inhomogeneous noise, in which case the inhomogeneity in 
the noise couples the different Fourier modes of x, thus 
leading to a problem in the field of partial differential 
equations, was studied in [23 |. 



